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We consider a binary mixture of colloid and polymer particles with positions on a simple cubic 
lattice. Colloids exclude both colloids and polymers from nearest neighbor sites. Polymers are 
treated as effective particles that are mutually non-interacting, but exclude colloids from neighboring 
sites; this is a discrete version of the (continuum) Asakura-Oosawa-Vrij model. Two alternative 
density functionals are proposed and compared in detail. The first is based on multi-occupancy 
in the zero-dimensional limit of the bare model, analogous to the corresponding continuum theory 
that reproduces the bulk fluid free energy of free volume theory. The second is based on mapping 
the polymers onto a multicomponent mixture of polymer clusters that are shown to behave as hard 
cores; the corresponding (Highlander) property of the extended model in strong confinement permits 
direct treatment with lattice fundamental measure theory. Both theories predict the same topology 
for the phase diagram with a continuous fluid-fee freezing transition at low polymer fugacity and, 
upon crossing a tricritical point, a first-order freezing transition for high polymer fugacities with 
rapidly broadening density jump. 

PACS numbers: 61.20.Gy, 82.70.Dd, 64.75. +g 



I. INTRODUCTION 

Mixtures of colloidal particles and non-adsorbing poly- 
mers suspended in a common solvent P, have been 
theoretically investigated on various levels of descrip- 
tion, ranging from simplistic to realistic effective inter- 
actions between the constituent particles. The proto- 
type of the former is the Asakura-Oosawa-Vrij (AOV) 
model 0, H that describes the colloids as hard spheres 
and the polymers as ideal (i.e., non-interacting) effective 
spheres that interact via hard core repulsion with the 
colloids. Despite its simplicity this model reproduces the 
essential trends in bulk phase behavior of colloid-polymer 
mixtures, involving colloidal gas, liquid, and crystalline 
phases 0, IE S 13 j an d proved to be useful for studying 
a range of interfacial properties 0, like the structure of 
colloidal gas- liquid interfaces and wetting of substrates, 
issues that are experimentally relevant 0, 0, . 

Density functional theory [liiL IT^H is a primary tool to 
treat spatially inhomogeneous systems. For the common 
reference model of additive mixtures of hard spheres, the 
fundamental measures theory (FMT) 15] is many inves- 
tigators' current choice. This approach was extended to 
cover the AOV model 0], thereby triggering much fur- 
ther interest in the study of interfacial properties of this 
model, see e.g. 0| f° r recent work on novel effects in 
sedimentation-diffusion equilibrium. In a different direc- 
tion, FMT was recently generalized to hard core lattice 
models 0, 0, |2(j elucidating the very foundation of the 
FMT approach. 



Despite its usefulness for practical applications, the 
AOV functional of Ref. |n| possesses several deficiencies, 
that are absent in the hard sphere FMT (see Ref. [2l| for 
a detailed discussion): i) When applied to one spatial 
dimension a spurious phase transition is predicted |2lj |. 
that is absent (as befits a model with short-ranged in- 
teractions) in the exact solution [22j . ii) In three dimen- 
sions the location of the bulk fluid-fluid critical point lies 
at lower polymer fugacity as compared to the value ob- 
tained with simulations, iii) The excess free energy is a 
linear functional of the polymer density profile. 

While we will not present a full resolution of the above 
issues, dealing with a lattice model allows us to go sig- 
nificant steps beyond the recipe of construction of Ref. 

Hence we consider a simplified version of the AOV 
model by constraining the position coordinates to an un- 
derlying lattice, for convenience taken to be of simple 
cubic symmetry. We are inspired by the fact that similar 
lattice models have been proven useful in soft matter re- 
search to make qualitative predictions, like e.g. for demix- 
ing j^, |3 and three-phase equilibria [25| . Also much 
vital attention has been recently paid to lattice models 
with short-ranged attractive interactions to address ad- 
sorption in disordered porous media [2|| IH, |3(| • 

We consider the case of small particles sizes, namely 
such that particles only repel other particles from near- 
est neighbor sites. As the polymers are ideal, each 
site may be occupied by more than one of those par- 
ticles. We compare two different DFTs for the binary 
AOV lattice model, both based on the lattice funda- 
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FIG. 1: (a) Illustration of the lattice AOV model realized in 
d = 2 spatial dimensions. Colloidal particles (dark rhombi) 
exclude their site and their nearest neighbor sites to colloids 
and polymers. Polymer particles (light rhombi) exclude col- 
loids from their site and their nearest neighbor sites, but 
can overlap freely with other polymers. Integers indicate 
the number of (overlapping) polymers at occupied sites, (b) 
Equivalent non-additive multi-component hard core mixture. 
Clusters of overlapping polymers in (a) are viewed as quasi- 
particles (light rhombi) that exclude their site to other poly- 
mer clusters, and their site and their nearest neighbor sites to 
colloids. Integers enumerate the different species of polymer 
clusters and correspond to the number of overlapping poly- 
mers on the lattice sites in (a). Colloids (dark rhombi) behave 
as in (a). 



dimensional simple cubic lattice 1 d . The interaction be- 
tween colloids is that of site-exclusion and nearest neigh- 
bor exclusion, corresponding to the pair interaction po- 
tential 



Vcc(r) 



oo if |r| < 1, 
otherwise, 



(1) 



where r 6 % d is the center-center distance between the 
particles. Colloids and polymers interact similarly: 



Vcp(T) 



oo if |r| < 1, 
otherwise, 



while polymers are ideal (non- interacting) , 
V pp (r) = 0. 



(2) 



(3) 



In essence this is a discretized AOV model with equal- 
sized components; see Fig. ^ f° r a sketch. 



III. DENSITY FUNCTIONAL THEORIES 



A. Overview of LFMT 



mental measure theory (LFMT) [H E HJ. The first 
approach allows for multi-occupancy of polymers in the 
zero-dimensional limit, and hence constitutes the lattice 
version of the (continuum) theory of Ref. 16]. The sec- 
ond approach relies on the introduction of polymer clus- 
ters as quasi-species that feature hard-core interactions 
only. We hence arrive at an extended model that pos- 
sesses a "Highlander" [^2 property: In an appropriately 
small cavity, there can only be one particle. Exploiting 
further the fact that LFMT can cope directly with small 
non-additivities [l^, we arrive at a theory for the AOV 
model that is exact when applied to one dimension; hence 
it does correctly predict the absence of phase transitions. 
We find that although the derivations of the two DFTs, 
as well as their appearance at first glance differ markedly, 
the results for bulk phase behavior in three dimensions 
are very similar, demonstrating the internal consistency 
of FMT. 

The paper is organized as follows. In Sec. [n] we define 
the lattice AOV model. Sec. 11111 is devoted to the con- 
struction of both density functional theories. In Sec. IIVI 
results for the bulk phase diagram are presented and we 
conclude in Sec. 



II. THE MODEL 

We consider a binary mixture of particles representing 
colloids (species c) and polymers (species p) on the d- 



As is customary in density functional theory [l3L Il4| , 
we express the Hclmholtz free energy functional as F = 
F[d + F cx , where the free energy of the binary ideal lattice 
gas is 

Fu\poPp] = E I>( r )KPiW) - 1], (4) 

i—c^p r££ 

with pi(r) being the occupancy probability of site r by 
particles of species i — c,p, and £ denoting the lattice, 
here C = If 1 . For the excess contribution to the total 
free energy, _F CX , we will in the following use two differ- 
ent implementations of LFMT. This theory permits to 
obtain an approximation to F ox for a lattice model from 
the exact solution obtained in finite (and small) sets of 
lattice sites, alternatively called nodes 0,0, E3|. The 
essential step determining the accuracy of the theory is 
the choice of maximal cavities |2f_| . Those are maximal in 
the sense that any further cavity taken into consideration 
can be obtained as an intersection of maximal cavities. 
The common choice (inherited from the continuum ver- 
sion of the theory |32j, l33l 1.3 11 ] ) for these maximal cavities 
is to set them equal to zero-dimensional (Od) cavities, 
as detailed below. Once chosen, the cavities uniquely 
determine the form of the functional under the require- 
ment that it yields the exact result when evaluated at Od 
density profiles pfj| | . For a general binary mixture (with 
species labeled by c and p) the final result possesses the 
form 

F C x[pc,Pp}= E a(C)F c c x [p c ,p p ], (5) 

C cavities 
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FIG. 2: Maximal cavities along the X axis (the ones along 
the Y axis are just 90° rotations of these ones) for the multi- 
occupancy (a) and the Highlander (b) DFTs. Dark nodes can 
be occupied either by colloidal or by polymer particles; light 
nodes can only be occupied by colloidal particles. 



where the summation runs over the set of maximal 
cavities and their nonempty intersections; the a(C) are 
uniquely determined integer coefficients [2(|, and F^. is 
the exact excess free energy functional of the system when 
confined to cavity C. 



B. Multi-occupancy version of DFT 

The definition of Od cavities for lattice models with 
only hard-core interactions is unambiguous: Any set of 
nodes that can accommodate at most one particle consti- 
tutes a Od cavity. For the present model, however, this is 
problematic: As polymers are ideal, even a single lattice 
site can be multiply occupied by polymers, hence no Od 
cavity can contain polymers, and the resulting functional 
would account for them only through the ideal part, en- 
tirely ignoring the colloid-polymer interaction. One way 
out is to relax the definition: Any set of nodes in which 
any two particles present necessarily overlap constitutes 
a Od cavity Hence, as opposed to the first definition, 
the presence of more than one particle is not excluded a 
priori. (For mixtures with only hard core interactions, 
both definitions are equivalent.) 

For the present model it is straightforward to show 
that according to the modified (multi-occupancy) defini- 
tion, any pair of adjacent nodes on the simple cubic lat- 
tice constitutes a maximal Od cavity (Fig. [3^): Such an 
arrangement can either hold one colloid or an arbitrary 
number of (overlapping) polymers. Adding an additional 
node destroys the Od character as (at least) two nonover- 
lapping particles will fit. Implementing the scheme of 
Ref. H3 with this choice of basic cavities we obtain 



(3F cx [p c , Pp ] =^$(r), 



(0) 



with j3 = 1/(&bT), where T is temperature and ks is the 
Boltzmann constant, and the (scaled) free energy density 



$(r) = Y,^° V (4 a) (r),n^(r)) 

_a=l 

-(2d-l)^ ov (p c (r) )ft) (r)), 
depending on weighted densities 

n[ a \r) = pi(r) + p;(r + e ), 



(7) 



c,p, 



(8) 



where e a denotes the unit Cartesian vector along the 
a = 1 ..... ti axis and 



aAOV 



(Vc ri P ) = (1 - jfc - 77 P ) ln(l - r? c ) + i] c (9) 



is the excess free energy of a Od cavity for the AOV 
model [lr|: we give a derivation of 4>q° V (Vc Vp) m the 
Appendix. Since <fio° V (VcjVp) depends linearly on r\ v as 



where 

<h(v) = V + ( 1 -V) In(l-ry), 
(|7|) admits the rewriting 

d 



(10) 



(11) 



(2d-l)0o(pc) 



n: =1 (i-^ a) )(i 



Ppln 



(12) 



(1 - Pc) 



2d-l 



where, for the sake of notational simplicity, we omit here 
and in the following the dependence of p c ,p p ,nc on 
r and have introduced the short-hand notation ni"^ = 
n\ '(r-e Q ). 

The linear dependence of the functional on p p (r) yields 
a particularly simple form of the Euler-Lagrange equa- 
tion for the polymer density distribution, 



Pp 



nti(i 



» 



)(i 



» 



(i 



\2d-l 



(13) 



being explicit in p p (i.e., independent on the right 
hand side of p p ). Given that we have eliminated 
p p (r) as a functional of p c (r), p c (r ± e Q ), a — 
1, . . . , d, and z p (r), the Legendre-transformcd functional 
0T cx [p c ,z p ] = E r $ eff (r), with 



$ off = p p ln(p p /z p ) - p p + 



(14) 



is an excess free energy functional for the effective one- 
component fluid of colloidal particles interacting through 
polymer depletion; p p (r) is obtained from (|13fl . The 
strength of the depletion interaction is controlled by 
z p (r), which in turn can be determined by an external po- 
tential acting on the polymers and by the polymer fugac- 
ity (A more detailed account of the procedure to obtain 
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functional for effective fluids can be found in Ref. l35h 
With the appropriate substitutions, we obtain 



$ eff = 



(2d-l)Mpc) 



nLi(i-^ a) )(i-"^) 



(1 - Pc) 



2d-l 



(15) 



This completes the prescription for the lattice analog of 
the continuum AOV functional of Ref. 0] . 



species, n = 1. Second, we will handle the full (infinite) 
number of quasi-particle species. Subsequently mapping 
the result back yields a DFT for the lattice AOV model. 



2. Binary mixture with one quasi-species 

For the binary mixture of colloids and 1-polymer clus- 
ters, with density fields P c { r ) and pi(r), respectively, ac- 
cording to Refs. 0, |2fl the excess free energy density 
<&(r) is given by 



C. Highlander version of DFT 

1. Strategy of polymer clusters as quasi-particles 

Here we elude the presence of several particles inside a 
cavity and aim at sticking to the stronger (Highlander) 
definition: There can be only one |3lj| (here: particle in 
the cavity). For this purpose, polymer ideality consti- 
tutes a major problem, as in the AOV lattice model each 
site occupied by a polymer can be occupied simultane- 
ously by an arbitrary number of further polymers. To cir- 
cumvent this problem we will hence map the model onto 
an extended one: We refer to a piling of n polymers at the 
same node as a "polymer cluster" and treat polymer clus- 
ters as quasi-particles of a set of new species labeled by n. 
Thus there is a one-to-one correspondence between each 
configuration of polymers of the original model and a cor- 
responding (unique) configuration of polymer clusters of 
the new model. (Colloids are treated as before.) We re- 
gard polymer clusters as having shapes smaller than the 
original polymer, see Fig.^D. As a consequence, polymer 
clusters behave as hard core particles with i) site exclu- 
sion to other polymer clusters and ii) site exclusion and 
nearest neighbor exclusion to colloids. Clearly, feature ii) 
is directly inherited from the colloid-polymer interaction. 
Feature i) may be unexpected at first glance. Consider 
that a site occupied by an n-polymer cluster cannot be 
simultaneously occupied by an m-polymer cluster. Al- 
though in the original model the site can well be occu- 
pied by to + n polymers, in the extended model such a 
configuration corresponds to one single (to + n)-polymer 
cluster. Hence polymer clusters repel like hard cores do. 

Despite being nonadditive, the new model belongs to a 
class of nonadditive lattice models which are amenable to 
LFMT and whose one-dimensional version was shown to 
be exact 0] . The Od cavities for the new model are non- 
trivial and consist of two adjacent nodes available for the 
colloids, one of which (but only one) is also available to a 
quasi-particle (polymer cluster) of any species (Fig. |5Jd) . 
Clearly, the presence of either one colloidal particle or 
one quasi-particle excludes other particles from the cav- 
ity, and such cavities are obviously maximal. 

Here we will construct the corresponding free energy 
functional in two steps. First, for simplicity, we con- 
sider the quasi-particle model with a single quasi-particle 



$(r) = ^2{^(4 a Hr)+ Pl (r)) 

i=l 

+0 O (nW(r - e a ) + pi(r)) - O (n^\v)) } 
-(2d-l)0o(/>c(r)+pi(r)), (16) 

with 0q (77) defined in Eq. (|llfl . Thus, given a local fu- 
gacity field zi(r), the Euler-Lagrange equation for the 
1-polymer cluster density is 



91 (I-Pc-P!)-- 1 ' (17) 

where again the dependence on r has been omitted. Eq. 
(|17fl is an implicit algebraic equation for pi(r), so it de- 
termines pi(r) as a functional of p c (r) and zi(r). This 
permits to obtain, as in the previous case, a functional 
/3TCox[pc, zi] for the effective one-component fluid. 



3. Mixture with an infinite number of quasi-species 

Incorporating the infinite number of quasi-species 
amounts to first correctly modifying the ideal free energy 
functional to (3F id [p c , {p n }] = J2 r *id(r), with 



$ w (r) =p c (r) ln(p c (r)) - p c (r) 

OO 

+ E [^n( r ) ln ( V nP"( r ) 



Pn(r) 



(18) 



where p n (r) denotes the density profile of the quasi- 
species of n-polymer clusters, and V„ a thermal "vol- 
ume" accounting for the internal partition function of 
the polymers inside the cluster. (It is not hard to guess 
that V n — n\; nevertheless, this result will emerge ex- 
plicitly from the subsequent analysis.) Secondly, we use 
again LFMT for the excess part of the free energy. The 
necessary modification over the binary case fSec. HII C 2() 
in order to include the infinite number of quasi-species is 
actually very simple : we replace pi (r) in <|16|l by the 
total density of polymer clusters (pc) , 



X 

E 



Pn(r). 



(19) 
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As our aim is to describe the lattice AOV model, we 
have to "project" the quasi-species back to polymers. 
Clearly, the total cluster density and the polymer den- 
sity are related via 



Pp 



(r) = J2np n (r) 



(20) 



For convenience we will again resort to an effective one- 
component description. Hence we transform to the semi- 
grand potential 



PT[p c , z p ] = mm (*id(r) + *(r) - p P (r) lnz p (r)j , 

(21) 

where $id(r) is given by (fl£)) . and $(r) by l(TH|l with pi(r) 
replaced by p pc (r). With the same notation as in (|17|l . 
performing the minimization in l|21|l yields 

_ z p nLiC 1 - n c a) - p P c)(i - - p P c) , > 
Pn "v„ (i-pc-ppc) 2 ^ 1 ' [ ' 

which using l|19|) leads to 



Ppc = C(z P ) 



nLi(i 



» 



Ppc)(l 



» 



Ppc) 



(1 - Pc - Ppc) 



2d-l 



(23) 



with C{z p ) = J2n=i( z p/ V n)- Also, from Eqs. J20I), 123 
and (f?5f> we get 



P P 



7 qM 

P C(Zp) 



Ppc 



(24) 



establishing a simple proportionality between p p (r) and 
Ppc(r)- 

The function C,(z p ) can be easily obtained realizing that 
it is independent of all densities appearing in ll'Mli , rather 
it depends solely on z p (r). Thus if we particularize (|2*3)l . 
e.g. for p c ( r ) = 0, we obtain the simple relationship 



Ppc = C{Zp){l - Ppc), 



from which 



Ppc 



C(z P ) 



(25) 
(26) 



Then, according to Eq. (|2H) . 

P, = ^T^)=%(M1 + C))' (27) 

and since in the absence of colloidal particles, the poly- 
mers form an ideal gas (hence p p = z p ), we get 



ln(l+C)J =1- (28) 
The solution to this equation satisfying ("(0) = is 



c(%) = ^-i = £tt> 



(29) 



providing the confirmation of the relationship V„ = n\. 



4- Functional for the lattice AOV model 

Returning to Eq. I|23[l . we observe that this equation 
for Ppc( r ) is formally equivalent to Eq. ((T7|l for pi(r) with 
zi(r) replaced by the uz p (r)) given by J2HJ). Recalling 
the origin of Eq. (|17|l , this means that we can describe the 
original AOV model in terms of the free energy functional 
PF[p c , Ppc ] = E r {<i>id(r) + $(r)}, with 

$id = PcQnpc - 1) +Ppc(hl/9p c - 1) 
d 

$ = [M4 a) + Ppc) + M4 a) + ppc) 

i=l 

- M4 a) )] - (2d - l)MPc + Ppc), (30) 

with the fugacity for the total density of polymer clusters, 
Ppc(r), given by 



C(r) 



o z p( r ) 



1. 



(31) 



and the polymer density obtained from p pc (r) as 

P P (r)= ^ M Ppc(r). (32) 



This completes the prescription of the Highlander func- 
tional for the AOV model. 

As a final note, we can alternatively obtain the effective 
functional f3T ex [p c , z p ] = J2 r ^ eS ( r )> wmcn , after a few 
algebraic manipulations, can be written as 

$ off = p c In p c - (2d - 1)(1 - p c ) ln(l - p c - Pp C ) 

d 

+ E[( 1 -o in (i-<-p P c) (33) 

Of — 1 ^ ' 

+ (l-n?)ln(l-n?-ft»c) 
-(l-<)ln(l-<)], 

with Ppc(r) being the solution of l|23|) . 

D. Relationship between both DFTs 

Assuming small polymer fugacity, z p (r) <C 1, for the 
Highlander DFT, from Eqs. ^ and it follows that 
C(r) ps z p (y) and p pc (r) ~ P P ( r )- The functional is then 
approximately given by Eq. i|16|) . On the other hand, 
small fugacities imply small densities, so we can expand 
(|16fl to linear order in p p (v). Doing so and taking into 
account that the sum over r allows us to shift the argu- 
ments of the functions, we realize that the expansion is 
precisely functional Q). So both DFTs coincide at low 
polymer fugacities. 



IV. RESULTS 

As an application we consider the bulk phase diagram 
for the lattice AOV model and compare the predictions 
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from both DFTs. Guided by its continuum counterpart, 
we expect the phase behavior of the lattice model to be 
determined by the interplay of hard core colloid-colloid 
repulsion and the (short-ranged) colloid-colloid attrac- 
tion induced by polymer depletion. Thus a clear frame- 
work to study the model is the effective one-component 
description, cf. Eqs. I|15|) and l|33|l . As in other (one- 
component) models with hard core repulsion and short- 
range attraction, both condensation (corresponding to 
demixing from the viewpoint of the binary mixture) and 
freezing may occur. We can study the former via a con- 
vexity analysis of the effective (semi-grand) free energy 
of a uniform fluid, and the latter by considering spa- 
tially inhomogeneous density distributions characteris- 
tic for crystalline states. Here the candidate is a face- 
centered cubic (fee) structure that constitutes the close 
packed state for the colloids. In the following we will 
carry out this program for either DFT. 



nodes of sublattices a and b are denoted by p a and pb, 
respectively. Hence the average colloid packing fraction 



Vc = Pa + Pb- 

The semi-grand potential is given by 

y = y lnp a + — bxpb + d(l - 77 c ) ln(l - rj c ) 

~(d-±) {1- Pa ) ln{l- Pa ) 

~{d-i) (l- Pb )Hl- Pb ) 

(1 - Vc ) 2d (1 - 77c) 2d 



(37) 



(38) 



{I- Pa) 



2d-l 



(1 - Pb) 



2d-l 



The equilibrium density distribution for given values of 
r) c can be obtained by minimizing Y with respect to p a 
and pb under the constraint l|37[) . leading to the condition 



Thermodynamics from the multi-occupancy 
DFT 



As an appropriate thermodynamic potential we con- 
sider a (scaled) semi-grand free energy per volume, Y = 
p c [ln(p c ) — 1] + /3T ex /V\ For a uniform fluid at given 
polymer fugacity z p , according to (|15|l this is given by 



F = |ln(|)+d(l-r ?c )Ml-r ?c ) 
_ 1 

c 

(1 - Vc) 2d 



( M -i)(i-f)m(i-f 



(34) 



where r\ c = 2/? c is the colloid packing fraction. Convexity 
implies Y"(r\ c ) > 0, so the spinodal condition, Y"{r) c ) = 
0, yields 



d(2d-l)r )c (l-r ]c ) 



\Vc) 



2d— 1 



(35) 



where we have used the polymer reservoir packing frac- 
tion, rf p = 2z p , as an alternative thermodynamic variable. 
The minimum of this curve determines the critical point, 
given by one of the roots of the cubic polynomial 

(d - l)(?7 c crit ) 3 + 2(d - l) 2 (r? c crit ) 2 + (2d + l)?7 c crit -2 = 0, 

(36) 

and the corresponding value ?7p' crlt obtained from (|35|l . 

For d = 3 the critical point is at 77™* = ^/5/2- 1 » 0.225 
(see Fig OJl. The liquid-gas binodal can be obtained from 
()34|l via a double tangent construction, in practice carried 
out numerically. 

For crystals we are guided by the close-packed state 
of the colloids, and distinguish two fee sublattices of the 
underlying simple cubic lattice, each formed by the near- 
est neighbor nodes of the other. The colloid density at 



Pa(l - Pa) 



2d-l 



exp 



-(2d- 1)^1 - 



Pa 



2d\ 



p b (l - PbY 4 - 1 exp -(2d-l)z„ 



Pb 



(39) 



This equation possesses one solution characteristic for 
fluid states, p a = Pb = Vc/2, which is locally stable as 
long as it minimizes Y(p a ,pb), he. up to the value of r) c 
at which 



d 2 

jr^y{p a ,Vc - Pa, 
dpi 



= 0, 



(40) 



defining the freezing spinodal, that we can explicitly ob- 
tain as 



2(l-f) rf (l-dr? e ) 
d(2d - 1)^(1 - r) c ) 2d 



(41) 



If freezing is second order, Eq. I|41|l directly gives the 
transition line in the phase diagram. In the case of a 
first-order transition again a numerical double tangent 
construction is required to obtain the coexistence densi- 
ties. 



B. Thermodynamics from the Highlander DFT 

In the fluid phase we obtain the semi-grand potential 

(D-px-^-fM 1 - !-«><) 

(1 - r) c - p pc )^j 

(42) 



Y =^ln 
2 



d(l - ?7 C ) In 



1 - Vc 
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with p pc being the solution of 

l-i 

1 - - Pnr ) Ppc = C(z p )(l - rj c - p pc ) 2 . (43) 



Vc 

2 



Again, the fluid-fluid spinodal can be determined by 
Y"(t) c ) — 0; necessary expressions for p' pc {r]c) and p pc {vc) 
can be obtained implicitly from l|43[) . The spinodal poly- 
mer cluster density is 



pc 



l (Vc) = 



(1 - 77c) [l + (d - l)/7c 



l + (d- l)(2d+l)»? c ' 
which inserted into (|43|) gives the liquid-gas spinodal 



(44) 



\ti c ) =2 In 1 



2 (2d - 1 



[4d(d-l)f 
[2d - 1 - (d - l)r? c ] 



[1 + (d - l)r/ c 



2d-l < 



%(1 - Vc) 



2d-l 



(45) 



[notice that rf p = 2 ln(l + £)]. The critical point is at the 
minimum value of 77p ,spin versus rj c , which we can obtain 
analytically as 



it l-d(d+l) + dV2-6d + 5d 2 
/c (d-l) 2 (2d+l) ' 1 ' 

In order to treat crystalline states we use the same 
density parameterization as in the treatment with the 
multi-occupancy version, with p a and pb being the densi- 
ties of two sublattices, and r/ c = p a + pb- The semi-grand 
potential is obtained as 

Y = Y lliPa + Y Pb 



-(d-^)(l-p a )\n(l-p a -p p a J) 
-(d-\)(l-p b )Ml-pb-p v h) c ) 



+ d(l -7? c )ln 



and ppc'', a = a,b, fulfill 

Pp i(i - Pa - Pp i ) 2d - 1 - c(*p)a - »?» - ? d . (48) 

Minimizing F(p a , ?7 C — Pa) with respect to p a leads to 

p a (l - p a - Pp a }f d - 1 = Pb {l -p b - $>) M - X . (49) 
We obtain dp p a c / dp a through implicit differentiation in 

(gHJi. 

As in section IIV Al the freezing spinodal is obtained 
through Eq. I|40|) . yielding 



1 - Vc 



(47) 



rpc 



Pa=Pb=Va/2 



(1 - Vc )(l - dr, c ) 
1 + d(2d-3)r] c 



(50) 



which combined with (|48|l gives the analytic result 

2(2d- l) 2 ^- 1 



2 In 



[4d(d-l)] 



2d 



(2d-l-d Vc ) 2d - 1 f 



(51) 



C. Phase behavior in one dimension 

It is clear from a general argument applicable to one- 
dimensional systems with short-ranged forces that the 
lattice AOV model does not exhibit a phase transition for 
d = 1. In this case the Highlander functional with 
p p (r) replaced by p pc (r), is exact ^8|. (Note that the 
mapping of the polymers to hard-core polymer clusters, 
to which we apply the theory of 0, is an exact trans- 
formation.) Hence an immediate consequence is that 
the Highlander version correctly predicts the absence of 
phase transitions. We can recover this result explicitly: 
the position of the liquid-gas critical point, obtained via 
is 77™* = 1/2, rjp — oo, consistent with the absence 
of the transition. The freezing spinodal, Q5ip. reduces to 
rf p — oo, again reflecting the absence of a phase transi- 
tions. The multi-occupancy version, however, incorrectly 
predicts both condensation and freezing, see Eqs. (|35|) 
and JHJ. The one fluid phase is below (in rfp ) the liquid- 
gas critical point located at rj c J lt — 2/3, ?7p' crlt = 4. 



D. Phase behavior in three dimensions 

We display the result for the full phase diagram from 
either DFT in Fig. [3] plotted as a function of colloid 
packing fraction, r/ c , and polymer reservoir packing frac- 
tion, ri p \ For rf v = the system is a pure (colloid) hard 
core lattice gas with nearest-neighbor exclusion. Both 
DFTs reduce to the (same) LFMT that predicts a con- 
tinuous freezing transition at a (colloid) packing fraction 
of 77 c = 1/3. This is considerably lower than the value 
obtained from Pade approximants [3(j, 77™' w 0.43. 

Increasing 77^ leads to a shift of the transition to smaller 
values of i] Cl which can be attributed to polymers sub- 
stituting colloids on crystal lattice positions and hence 
decreasing the colloid density at the transition. At a 
threshold value of 77^ the transition becomes first order 
and a density gap opens up. The location of this tri- 
critical point differs somewhat in both treatments, being 
located at a larger values of rf v in the Highlander DFT. 
Upon increasing 77^ further, the coexistence density gap 
increases strongly. The liquid-gas transition is found to 
be metastable with respect to freezing. (The liquid-gas 
binodal obtained from the Highlander version is again 
located at higher values of 77^ as compared to the multi- 
occupancy version.) 

The occurrence of a tricritical point for the freezing 
transition is a peculiarity of the lattice model absent in 
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FIG. 3: Bulk phase diagram of the lattice AOV model as a 
function of colloid packing fraction, r/ c , and polymer reservoir 
packing fraction, rf p . (a) Result from the multi-occupancy 
DFT. For small rf p , below the tricritical point (triangle), there 
is a continuous freezing transition into an fee crystal; its loca- 
tion for rf p = is at a smaller value of r/ c as compared to the re- 
sult by Gaunt [3(| (arrow). Above the tricritical point freezing 
is discontinuous; coexistence is along horizontal tielines (dot- 
ted lines) between the fluid and crystal branches of the bin- 
odal (solid lines). Also shown is the freezing spinodal (long- 
dashed line) where the fluid phase loses its metastability upon 
increasing r/ c . The gas-liquid binodal (short-dashed line) end- 
ing in a lower critical point (filled dot) is metastable with re- 
spect to freezing, (b) The same as (a), but obtained from the 
Highlander DFT. For comparison the tricritical point (open 
triangle) and the metastable gas-liquid critical point (open 
dot) obtained from the multi-occupancy DFT are shown. 



proaches the hard sphere crystal-fluid coexistence den- 
sities in the continuum model (without an intervening 
tricritical point). 

Given the fact that we do find the correct (hence 
also experimentally observed) fee crystalline structure, 
we hence conclude that the current model is very suit- 
able for the study of inhomogeneous situations at high 
rf pl where a dilute colloidal gas and a dense crystal are 
the relevant bulk states. 



V. CONCLUSIONS 

We have carried out a detailed comparison between 
two density functional theories for a lattice AOV model 
consisting of a binary mixture of colloidal particles and 
polymers both with position coordinates on a simple cu- 
bic lattice. The (pair) interactions are such that colloids 
exclude their site and nearest neighbors to both colloids 
and polymers, and polymers do the same with colloids, 
but do not interact with other polymers. Relying on 
the LFMT concept we have obtained two density func- 
tionals for arbitrary space dimension d by starting with 
the zero-dimensional Statistical Mechanics of the model. 
The multi-occupancy version is an analog of the DFT 
for the continuum AOV model . The Highlander ver- 
sion exploits the possibility of mapping the lattice model 
to a cluster model that features hard core interaction 
only. Both DFTs are exact for d = by construction. 
The Highlander version is also exact in d = 1, where the 
multi-occupancy version incorrectly gives a phase transi- 
tion. 

The predictions for the phase diagram for d — 3 from 
both approaches are very similar, and we have argued 
that the topology resembles that of real colloid-polymer 
mixtures for small polymer-to-colloid size ratios, where 
the liquid-gas transition is metastable with respect to 
freezing into an fee colloid crystal and the coexistence 
density gap is very wide for large values of polymer fu- 
gacity. We hence conclude that both the model and 
DFT approaches are well suited to study properties of 
inhomongeneous situations in colloid-polymer mixtures 
driven by dilute colloidal fluid and dense colloidal crystal 
phases. 



the continuum version, where freezing is first order for 
all rfp [?J. The phase behavior of the continuum AOV 
model depends sensitively on the polymer-to-colloid size 
ratio, q = a p /a Cl where a p and a c is the polymer and col- 
loid diameter, respectively. The phase diagram obtained 
for small values of q ~ 0.1 indeed roughly resembles the 
topology of the phase diagram that we find for the lattice 
AOV model, provided rL is sufficiently high. For small 
values of rf the broad coexistence region smoothly ap- 
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APPENDIX 

For notational simplicity, we use 1 = r and 2 = r + e Q . 
If Zi(x) denotes the local fugacity at node x(= 1,2) of 
species i(= c,p), then the grand partition function of a 
maximal Od cavity will be 



S = e *»(i)+Ma) + Zc (i) +Zc (2). 



(52) 



Since Pi(x) = Zi(x)<5 ln^o/<^(x), from the equation 
above 



Pc(x) 



z c (x) 



p p (x) = ^ Xje 



—0 



(53) 
(54) 



From J5 



- e>\ z c (l) + z c (2) 
Tic = Pc(l) +Pc(2) = , (55) 



so eliminating z c (l) + z c (2) in i|52|) we obtain 

e «,(i)H(>) = s (l-i ?c ). (56) 
With this Eq. ij5*4*)l becomes 

p p (x) = Z p (x)(l - T] c ), (57) 

what allows us to obtain from 156(1 



InSn = 



l-ric 



ln(l - r, c ) 



(with the obvious notation rj p = p p (l) + p p (2)). 
Now, as 4>q ov = [3F cx [p c , p p ] in this cavity, 



4 aov 



(58) 



In S + Pi( x ) ln (W x )/M x )) > ( 59 ) 
and using , J37|) and JSHJ we finally obtain © . 
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